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Abstract. We review recent work of the authors on the non-relativistic Schrodinger 
equation with a honeycomb lattice potential, V. In particular, we summarize results 
on (i) the existence of Dirac points, conical singularities in dispersion surfaces of Hy = 
— A + V and (ii) the two-dimensional Dirac equations, as a large, but finite time, effective 
description of e~ lHvt i/jo, for data ipo, which is spectrally localized at a Dirac point. We 
conclude with a formal derivation and discussion of the effective large time evolution 
for the nonlinear Schrodinger - Gross Pitaevskii equation for small amplitude initial 
conditions, ipo. The effective dynamics are governed by a nonlinear Dirac system. 

o 
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O 

q ! 1. Introduction 

There has been intense interest within the fundamental and applied physics communities 
in the propagation of waves in honeycomb structures. Two areas where such structures 
have been explored extensively are (i) condensed matter physics and (ii) photonics: 

IT) , 

(1) Graphene, a single atomic layer of carbon atoms, is a two-dimensional structure 
with carbon atoms located at the sites of a honeycomb structure with remarkable 
electronic properties; see, for example, the survey articles |13|. 114]. 

(2) Photonic / electro-magnetic propagation (linear and nonlinear) has been studied 
in dielectric honeycomb structures [91 [151 El E] • 

The basic mathematical model is a wave equation, defined on the 2-dimensional plane 
. with a medium whose material properties vary periodically according to a honeycomb 

pattern: Schrodinger's equation with honeycomb lattice potential, in the quantum set- 
ting, and Maxwell's equations, with honeycomb structured dielectric parameters, in the 
. electromagnetic setting. 

We focus on the time-evolution of the Schrodinger equation 
^. idtil> = (-A + V(x))il; = ifW, (1-1) 

where F(x) is a smooth, real-valued and periodic potential, defined on K 2 . Denote by 0, 
a choice of elementary period cell of V(x). We are interested in properties of the time- 
evolution, which relate directly to the special symmetry of the honeycomb lattice. To 
represent the time evolution operator, e~ lHvt , we quickly review the basic Floquet-Bloch 
spectral theory of periodic potentials. 

Consider the family of k— pseudo-periodic eigenvalue problems: 

H v $(x;k) =/i(k) $(x;k), x£l 2 , (1.2) 

$(x + v;k) = e ik ' v $(x;k), v G A. (1.3) 
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We introduce the family of k— pseudo-periodic L 2 functions, L 2 /. 

Li = {fGL 2 (n) : /(x + v) = e lk ' v /(x), v G a} . 

Here, 

(1) A denotes the period lattice for V(x). We denote by 

Ah = Zvi © Zv2 , where vi 

is the period lattice of the regular honeycomb structure; see Figure [TJ and 

(2) k— varies of the Brillouin zone, B, a choice of fundamental cell in the dual period 

lattice, A*. For honeycomb lattice potentials, 

/ l \ / i 

4"7T / 2 \ 4tt 2 



A^ = Zki © Zk2 , where ki ,_ . . . ,_ . 

V3 \ / V3 I _yj? 

2 / V 2 

and /3 can be chosen to be a regular hexagon, centered at the origin. The vertices 
of B fall into two equivalence class: points of K type and points of K' type, which 
alternate around the perimeter of £>; see Figure [2j 

An equivalent spectral problem with periodic boundary conditions may be obtained by 
setting <l?(x;k) = e* kx p(x;k). Then, p(x;k), k G B is periodic and satisfies boundary 
value problem: 

ffy(k)p(x;k) = /i(k)p(x;k), x G M 2 , (1.4) 
p(x + v;k) = p(x;k), v G A, (1.5) 

where 

H v (k) = - (V + ik) 2 + V(x) . (1.6) 
For each kGB, the spectrum of iJy(k) is discrete, consisting of real eigenvalues: 

Mi(k)<Mk) <---<M&(k) <••• 

(listed with multiplicities), with corresponding eigenfunctions p&(x; k), 6 > 1. As k varies 
over i3, //&(k) sweeps out a closed real interval , fJ-b(B). The spectrum of Hy is the union 
of these intervals: 

spectrum (Hy) = (j E b (B) (1.7) 
6>1 

The Floquet-Bloch modes 3>b(x; k) = e* k x pb(x; k), b > 1 are complete in the sense that 
for any / G L 2 (R 2 ) 



/(x) - E / <<M-;k),/(-)} i2(R2) $ 6 (x;k)dk -> o 



(1.8) 



Kb<N 

"2 /"IB 2 



in 2/(IR 2 ) as N f oo. The solution of the initial value problem for time-dependent 
Schrodinger equation, fjl . 1[) . with initial data V( x , 0) = V'o(x), has the representation 

e- iHt ^o = E / e_W(k) * (^(•;k),^o(-))L 2 (R 2 ) * 6 (x;k)dk. (1.9) 
6>1 Jb 

An understanding of the time-dynamics (jl.9p requires, in particular, a detailed under- 
standing of the functions /ife(k), 6 > 1. These are called band dispersion functions of Hy. 
The graphs k h-» //&(k) are called the dispersion surfaces of Hy. 
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Figure 1. Part of the honeycomb structure, H. H is the union of two 
sub-lattices Aa = A + Ah (blue) and Ab = B + A^ (green) . The lattice 
vectors {vi,V2} generate A^. 




-0.5 0.5 



Figure 2. Brillouin zone, B, and dual basis {ki,k2}. K and K' are la- 
beled. Other vertices of Bh obtained via application of R, rotation by 2tt/3. 



Let's now introduce the class of potentials that interests us. For any function /, defined 
on R 2 , introduce 

ft[/](x) = /(iTx), (1.10) 
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where R is the 2x2 rotation matrix, which clockwise-rotates a vector by an angle of 2ir/3. 

Definition 1.1 (Honeycomb lattice potentials). Let V be real-valued and V G C°°(M. 2 ). 
V is a honeycomb lattice potential if there exists xo £ M 2 such that V"(x) = V(x — Xo) has 
the following properties: 

(1) V is A h - periodic, i.e. V(x + v) = V(x) for all x£l 2 and v € A^. 

(2) V is even or inversion-symmetric, i.e. V{— x) = V"(x). 

(3) V is 1Z- invariant, i.e. 

rc[v](x) = y(i?*x) = y( x ), 

where, i?* is the counter-clockwise rotation matrix by 2ir/3, i.e. R* = . 

Thus, a honeycomb lattice potential is smooth, A^- periodic and, with respect to some 
origin of coordinates, both inversion symmetric and 1Z- invariant. 

In this article we summarize recent results of the authors on the spectral properties 
of Hy and properties of the evolution operator, e ~ lH vt^ where V denotes a honeycomb 
lattice potential. In section [2] we describe results on the existence of conical singularities 
in the dispersion surfaces of honeycomb lattice potentials for quasi-momenta located at 
the vertices of B. In section [3] we consider the initial value problem for the non-relativistic 
Schrodinger equation, (jl.ip . with a honeycomb lattice potential. We study the evolution 
for initial conditions which are spectrally localized about a vertex (K- type or K'— type 
point) of B and show that the large, but finite, time-evolution is effectively governed by 
the constant coefficient two-dimensional Dirac equation. Finally in section 2] we remark 
on some open problems concerning the relation of the nonlinear Schrodinger / Gross- 
Pitaevskii equation and the two-dimensional nonlinear Dirac equation. 

Acknowledgements: The authors wish to thank M. Ablowitz, A.C. Newell and G. 
Uhlmann for stimulating discussions. 

2. Dirac points 

In this section we discuss results of the authors in [7] on the conical singularities, so-called 
Dirac points, in the dispersion surfaces of Hy, where V is a honeycomb lattice potential. 

A key property of honeycomb lattice potentials, V is that if K* denotes any vertex of 
B, then we have the commutation relation: 

[K,H V (K*)} = . (2.1) 

Since 1Z has eigenvalues 1 , r and r , it is natural to the split , the space of K*- pseudo- 
periodic functions, into the direct sum: 

Lk = L k,i © L h,r e L 2 Ki!>r , (2.2) 

where L 2 ^ a are invariant eigen-subspaces of 1Z, i.e. for a = l,r, r, where r = exp(27rz/3), 
and 

= {se L 2 ^ : Tig = ag] . (2.3) 
We next give a precise definition of a Dirac point. 
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Definition 2.1. Let V(x) be a smooth, real- valued, even (inversion symmetric) and pe- 
riodic potential on M 2 . Denote by B, the Brillouin zone. We call K 6 B a Dirac point if 
the following holds: There exist an integer b\ > 1, a real number /i*, and strictly positive 
numbers, A and 5, such that: 

(1) /U* is a degenerate eigenvalue of H with K— pseudo-periodic boundary conditions. 

(2) dim Nullspace n*l\ = 2 

(3) Nullspace (^H - = span j<3?i(x), <3? 2 (x)j, where <3?i G r and $2(x) = 

(4) There exist Lipschitz functions /i-t(k), 

M6i(k) = M-(k) /i& 1+ i(k) = /u + (k) 5 //±(K) = ii± 

and -E'-i-(k), defined for |k — K| < S, and k— pseudo-periodic eigenfunctions of H: 
<&±(x;\i), with corresponding eigenvalues /U-t(k) such that 

M + (k)-/i(K) = +A |k-K| (1 + £7+(k) ) and 

M_(k)-/i(K) = - A |k-K| (1 + S_(k)), (2.4) 

where |.E±(k)| < C|k - K| for some C > 0. 



Remark 2.1. In [7] we prove the following 

Proposition 2.2. Suppose conditions l.,2. and 3. of Definition {2A\ hold and denote by 
{c(m)} mgl 5 the sequence of L^ T Fourier- coefficients o/ < I ) i(x). Define the sum 

A a = £ c(m) 2 C M • KT (2.5) 
//Ajj 7^ 0, then 4. of Definition ^. 1\ holds (see (|2.4p J w;zi/i A = |Ajj|. 

Therefore Dirac points are found by verifying conditions l.,2. and 3. of Definition 12.11 
and the additional (non-degeneracy) condition: Aj ^ 0. We use this characterization to 
prove the following result, Theorem 5.1 of [7], concerning the existence of Dirac points for 
Schrodinger operators with a generic honeycomb lattice potentials: 

Theorem 2.1. Let V(x) honeycomb lattice potential. Assume further that the Fourier 
coefficient of V, Vn, is non- vanishing, i.e. 

V lfl = [ e-^+^-y V(y) dy + . (2.6) 
Jn 

Consider the one-parameter family of honeycomb Schrodinger operators defined by: 

ffte = _A + e U(x) . (2.7) 

There exists a countable and closed set Cct such that for all e ^ C, the vertices, K*, of 
Bh are Dirac points in the sense of Definition 12.11 

More specifically, the following holds for e ^ C: There exists b\ > 1 such that fj^ = 
Hi = +1 (K*) is a K*— pseudo-periodic eigenvalue of multiplicity two where 
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(1) [i% is an T - eigenvalue of of multiplicity one, and corresponding eigen- 
function, $|(x). 

Hi is an L^- f - eigenvalue of H^ e > of multiplicity one, with corresponding eigen- 
function, $|( x ) = ^i( _x )- 

\i% is not an x - eigenvalue of . 

(2) There exist 5 e > 0, C e > and Floquet-Bloch eigenpairs: (<J>+(x; k), fi e + (k)) and 
($l(x; k), /zl(k)), and Lipschitz continuous functions, E±(k), defined for |k — 
K*| < 5 e , such that 

4(k)- M £ (K,) = + |AJj| |k-K*| (1 + E £ + (k)) and 

/4(k)-//(K*) = -|AJ| |k-K*| (l + £L(k)), 

where 

\\ = J>(m,//,e) 2 ( M - K- / (2.8) 

is given in terms of {c(m, //, e)}, the L 2 ^ - Fourier coefficients of <3? e (x; K*). 

Furthermore, |Ey_(k)| < C e |k — K*|. Thus, in a neighborhood of the point (k, fj,) = 
(K*,//£) E M 3 , the dispersion surface is closely approximated by a circular cone. 

(3) There exists e° > 0, such that for all e G (-e°, e°) \ {0} 

(i) eV\ i > ==> conical intersection of I s * and 2 nd dispersion surfaces 

(ii) eV\ \ < ==> conical intersection of 2 nd and 3 rd dispersion surfaces . 

Remark 2.2. In a forthcoming article, we present a general analytic perturbation the- 
ory of deformed honeycomb lattice Hamiltonians, for perturbations which commute with 
inversion composed with complex conjugation. Conical (Dirac) points persist for small 
perturbations of this type, although the conical singularities typically perturb away from 
the vertices of B. These results extend those of [7j and, in particular, include the case of 
a uniformly strained honeycomb structure. 

Remarks on the proof of Theorem \2.1\ [7J : Following the strategy outlined in Remark 12.11 
we first show, for all e sufficiently small and non-zero, by a Lyapunov- Schmidt reduction, 
that the degenerate, multiplicity three eigenvalue of the Laplacian with K*— pseudo- 
periodic boundary condition, splits into a multiplicity two eigenvalue and a multiplicity 
one eigenvalue, with associated 1Z- invariant eigenspaces. Next we must show the persis- 
tence of this degenerate multiplicity two subspace as e is increased without bound. To 
continue to e large we introduce a globally-defined analytic function, £(/j,,e), whose zeros, 
counting multiplicity, are the eigenvalues of H( e >, Eigenvalues occur where an operator 
I + C(fi, e), C(fi,e) compact, is singular. Since C((j,,e) is not trace-class but is Hilbert- 
Schmidt, we work with £(fj,,e) = det2(/ + C(fi, e)), a renormalized determinant. Next 
£(/j,,e) and Ajj are studied using techniques of complex function theory to establish the 
existence of Dirac points for arbitrary real values of e, except possibly for a countable 
closed subset, Cel. The origin of this exceptional set is a topological obstruction to the 
existence of continuous choice of eigenvector as parameters are varied. 

Previous analyses of such honeycomb lattice structures are based upon extreme limit 
models: the tight-binding / infinite contrast limit (see, for example, [18 } 113 } [TT]) in which 
the potential is taken to be concentrated at lattice points or edges of a graph; in this limit, 
the dispersion relation has an explicit analytical expression, or the weak-potential limit, 
treated by formal perturbation theory in [9J, [2] and rigorously in [8]. 
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3. Evolution of wave packets and the two-dimensional Dirac equation 

From the previous section, we see that generic honeycomb Hamiltonians, Hy, to have Dirac 
points, and that some band dispersion function is approximately that of a two-dimensional 
wave equation with wave speeds =t|A(j|. Thus, initial conditions which are spectrally local- 
ized about a Dirac point, should evolve approximately according to a constant coefficient 
hyperbolic PDE with wave speeds i|Aj|. In this section we summarize results deriving the 
appropriate wave equation and its time-scale of validity in terms of the degree of spectral 
localization, 5, about a Dirac point, K*. 

A general solution of the time-dependent Schrodinger equation, constrained to the de- 
generate 2-dimensional eigenspace associated with eigenvalue, = /z(K*), associated 
with the Dirac point, K*, is of the form 

V>(x, t) = e-^ ( a x $i(x) + a 2 $ 2 (x) ) , (3.1) 

where a\ and ct 2 are arbitrary constants. 

Consider now a wave packet initial condition, which is spectrally concentrated near K*: 

^ 5 (x,0) = ^o(x) = 6 ( a 10 (<5x) ^(x) + a 20 (<5x) $ 2 (x) ) 

= 5 ( a 10 (5x) Pl (x) + a 20 (5x) p 2 (x) ) e lK - x (3.2) 

Here, 5 is a small parameter. We assume aio(X) and a 2 o(X) are Schwartz functions of 
X. We expect that this assumption can be weakend considerably without difficulty. The 
overall factor of 5 in (|3.2|) is not essential (the problem is linear), but is inserted so that 
ipQ has L 2 (M?)- norm of order of magnitude one. 
We seek solutions of (|1.1|) . (|3.2|) in the form: 

^(x,t) = e" 1 ^ 1 6 aj(«5x,«Jt)$j(x) + n <5 (x,i)j . (3.3) 

where ^(x.O) = 0, ay(X,0) = a i0 (X), j = 1,2 (3.4) 
to ensure the initial condition (|3.2p . 

The goal is to show that the Schrodinger equation (jl.ip has a solution of the form (|3.3j) 
with an error term, r/ <5 (x, t), which satisfies 

sup || r] s (;t) ||^ (R2) = 0(5 T *), 5^0 . (3.5) 

0<t<p<5- 2+£ i 

for some > 0, provided the slowly varying amplitudes ay(<5x, 5t), j = 1,2 evolve ac- 
cording to the system of Dirac-type equations (j3.6|) - (|3.7j) . Here, p > and e± > are 
arbitrary. 

In [6] we prove 
Theorem 3.1. Assume 

5 °< x > - ( ) e i^i 2 

and let a(X, T) denote the global-in-time solution of the Dirac system 

cVai(X,T) = -A] (d Xl +id X2 ) « 2 (X,T) (3.6) 

d T a 2 (X,T) = -A„ (0^-10*,) «i(X,T) , (3.7) 

where ^ Ah £ C. loii/i initial data a(X, 0) = c?o(X). 

Consider the time- dependent Schrodinger equation, (jl.ip . where V(x) denotes a poten- 
tial for which the conclusions of Theorem \2.1\ hold, e.g. V(x) = eV/ l (x), where Vh is a 
honeycomb lattice potential satisfying V\ \ ^ and e is no£ in i/ie countable closed set C. 
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Assume initial conditions, ipo, of the form (|3.2|) , Fix any p > 0,£i > and s > 1. Then, 
/ias a unique solution of the form (|3.3|) . where for any \a\ < N 



sup 

0<t<p 5- 2+£ i 



»7 4 (- } t) / ^ = (<T*), 5^0 (3.8) 



/or some r* > 0. 



Remark 3.1. In connection with Remark 12. 2\ we also consider the analogous question 
of the dynamics of solutions for wave-packet initial data, spectrally concentrated at a 
Dirac point of the deformed honeycomb structure. In this case, the large, but finite, time 
dynamics are being give by tilted- Dirac equations. The latter can be mapped to the 
standard 2D Dirac equations by a Galilean change of variables. 

Remarks on the proof of Theorem \3.1[ We seek a solution of the initial value problem 
with wave-packet initial condition ()3.2p of the form (|3.3p . a slow space/time modulation 
of the degenerate subspace plus an error term r/ S (x,t). The Dirac equations ()3.6p - (|3,7p 
arise as a non-resonance condition, which ensures that rj S (x,t) is small on a time interval: 
< t < 0{5- 2+ ). Estimation of rf , via its DuHamel integral equation, requires a careful 
decomposition of the propagator, e~ l( -~ A+v ^ and analysis of its action on functions with 
quasi-momentum components supported near K*, a vertex of B, and those with quasi- 
momentum components supported away from K*. The resonant terms, which are removed 
by imposing equations ()3.6p - (|3.7p . arise from quasi-momenta near K*. A detailed expan- 
sion of the normalized Floquet-Bloch modes for such quasi-momenta is required. Such 
modes are discontinuous at K*. Components corresponding to quasi-momenta away from 
K* are controlled, via Poisson summation and integration by parts with respect to time, 
by making use of rapid phase oscillations in time. 

Formal derivations of Dirac-type dynamics for honeycomb lattice structures are dis- 
cussed in the physics [13] and applied mathematics [2] literature. A rigorous discussion of 
the tight-binding limit is presented in [1]. Conical singularities have long been known to 
occur in Maxwell equations with constant anisotropic dielectric tensor; see, for example, 
|10| . [1] and references cited therein. 



4. NLS / GP and the nonlinear Dirac equation 

A model of fundamental importance in the description of macroscopic quantum phenom- 
ena and nonlinear optical phenomena is the nonlinear Schrodinger / Gross-Pitaevskii 
equation, (NLS/GP): 

iM = (-A + y(x))^+ g\M 2 i> (4.1) 

We refer to the case V = as the nonlinear Schrodinger equation (NLS). In the quantum 
setting, the potential, V(x), models a magnetic trap confining a large number of quantum 
particles. A wide range of potentials, so called optical lattices, can be induced optically 
via interference of optical beams. We consider the case where V = Vh, a honeycomb 
lattice potential in the sense of Definition 11.11 Analogously, in the setting of nonlinear 
electromagnetics interference of beams in photorefractive crystals can generate a spatially 
dependent index of refraction, corresponding to a honeycomb lattice potential. See, for 
example, [El H2 E] ■ 

Our goal in this section is to give a formal derivation of a nonlinear Dirac equation from 
the NLS/GP and to pose some interesting questions to consider going forward. 
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As earlier, we assume wave-packet initial conditions for NLS/GP, with amplitude scaled 
so that nonlinear effects and the Dirac dynamics enter on the same time-scale: 

^(x,i) = ( 55j^a J -(X,T)$ i (x) + r/ 5 (x,i) j , (4.2) 

where X = 5x and T = St. Substitution of the Ansatz (g2]) into NLS/GP, (jH]), yields 
rf = 52 {5rji + 0(S 2 )), where r/i, satisfies: 



id t m (x, t) - (H - fj,+) r/i (x, t) 



= id T<Xj $i( x ) + 2V x aj • V x $j(x) - 5 2 "i^^ *i( x )*i( x )*K x ) 

3=1 \ Ki,l<2 

(4.3) 

A calculation which is analogous to that carried out and rigorously justified in [6] for 
the linear Schrodinger equation (see section [3]) yields the following semilinear Dirac system 
governing the slowly modulating amplitudes a%(X.,T) and a 2 (X,T): 

d T ai = -Ajj (d Xl +idx 2 ) a 2 - i g (/3i|«i| 2 + 2/3 2 |a 2 | 2 ) a\ , (4.4) 

9 t q;2 = -Ajj (d Xl - idx 2 ) cti - i g (2fa\a 1 \ 2 + Pi\a 2 \ 2 ) a 2 ■ (4.5) 

The coefficients j3\ and fa are given in terms of the Floquet-Bloch modes: 



fa 



( |$i(x)| 4 dx = f |$ 2 (x)| 4 dx and fa = I |$ 1 (x)| 2 |$ 2 (x)| 2 dx . (4.6) 
Jn Jn Jn 



The system ()4.4p - (|4.5p has the following conserved integrals: 
M[ai,a 2 ] = / ( I cti 1 2 + |ai| 2 ) dX\ dX 2 , 



% [«i, a 2 ] = Im 



Ah 



3 
4 



«2 (9xi +i0x 2 )«i dX 2 



fa |ai| 4 + 4/3 2 |ai| 2 |a 2 | 2 + fa |a 2 | 4 dXi dX 2 . 



(4.7) 



(4. 



A word about the calculation leading to (|4.4)) - (j4.5f) : Heuristically speaking, the mod- 
ulation equations (|4.4|) - (|4.5p for Oj(X,T), j = 1,2, are obtained by choosing their slow, 
(X, T), evolution to cancel all resonant components of the forcing term in (|4.3p . This can 
be implemented by setting the projection of the source term onto the modes $j, j = 1, 2 
(integrating only over the fast variable, x) to zero. 

In computing this projection, many cubic terms drop by symmetry considerations. 
For example, consider the contribution of the cubic term: aiaia 2 <&i&i& 2 in (|4.3p to the 
equation for a\. To compute this contribution take the L 2 (M 2 .) inner product with <3?i(x): 



($i,$i$i$ 2 ) 



L 2 (n) 



3>i(x) • $i(x)$x(x)$ 2 ( x ) dx 



_x=_R*y 



$i(fl*y) • $i(i?*y)$i(ii*y)$ 2 ( J R*y) dy 



R.n 



/ r*i(y) -r$i(y) -r$i(y) -T$ 2 (y) dy 

JRCl 

r 2 / *i(x) • $i(x)$i(x)l^(x) dx = r 2 ($i,$i$i$^) L2 



(O) 
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Therefore, 

(l-^K*i.*i*i^ ( n) =° <*i>*i^> i2( a)= - 
It follows that a term of the form aia^ai does not appear in the (|4.4p . 

An alternative formal approach to obtaining f)4.4|) - (|4.5j) is by the method of multiple 
scales; see [2] where a formal derivation of a semilinear Dirac system for a small ampli- 
tude ("shallow") honeycomb lattice potential of a special form is presented and studied 
numerically. 

As in the proof of Theorem 13. 1\ a proof of the large time validity of (|4.4|) - (|4.5p requires 
that i] S remains small in a suitable norm on a large (growing as 5 4- 0) time scale. This 
rests on a study of the time-evolution of the Floquet-Bloch components of rj s which are 
near K* and away from K*, together with an understanding of bounds on the solutions 
to the nonlinear Dirac equations (|4.4p - (|4.5|) . 

We conclude with some natural questions. Is there global well-posedness for the initial 
value problem for semi-linear Dirac system (|4.4p - (|4.5p or do solutions develop singularities 
(blow up) in finite time? How do the corresponding solutions of NLS / GP behave? 
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